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1. Introduction 

The Schrodinger functional (SF) is a useful tool to solve renormalization problems in lattice 
gauge theories [1^1]. A drawback of the standard formulation of the SF for Wilson fermions [Q] 
consists in large bulk 0(a) effects even in the massless limit. While these 0(a) effects can be can- 
celled by the usual Sheikholeslami-Wohlert term and 0(a) counterterms to the composite fields [|], 
§|, it is surprising that these bulk 0(a) counterterms are required at all, given the fact that Wilson 
fermions at zero mass enjoy the property of automatic 0(a) improvement [8|]. The origin of this 
problem is the explicit breaking of chiral symmetry by the standard SF boundary conditions for the 
fermions. Rendering automatic 0(a) improvement compatible with SF-like boundary conditions 
has been shown to be possible for theories with even flavour numbers [|| [J. An attractive solution 
are chirally rotated SF boundary conditions, which, for Nf = 2 flavours take the form 

Q+v(x)\x =o = = Q-y(x)\ Xo = T , y(x)QA =o = = y(x)Q-\x =t, (1-1) 

with Q± = ^(1 ±/yo/5T 3 ) and Pauli matrices t 1 - 2,3 acting in flavour space. The non-trivial flavour 
structure implies that 75 T 1 commutes with Q±, and can be used to resurrect the argument of auto- 
matic 0(a) improvement. Chirally rotated SF boundary conditions derive their name from the fact 
that they arise from the standard SF boundary conditions by performing a non-anomalous chiral 
rotation of the flavour doublet fields, 

Y -» R(oc)y, y \j/R(a) , R(a) = exp(/ay 5 T 3 /2). (1.2) 

The rotated fields satisfy boundary conditions involving the projectors P±(cc) = 4 [l ± 7oexp(/7sT 3 )] 
which interpolate between the standard SF boundary conditions (P±(0) = P±) and the chirally ro- 
tated ones in Eq. ( |1.1[ ), as P±(n/2) = Q±. Since chiral rotations are symmetries of the massless 
QCD bulk action, one may derive universality relations between correlation functions calculated at 
different values of a, 

(0[ ¥ ,xff]} m = (0[R(-n/2) V ,yR(-n/2)}) {P±) . (1.3) 

Here we have indexed the correlation functions by the projectors appearing in the boundary con- 
ditions. On the lattice with Wilson quarks and the standard Wilson bulk action, one expects to 
recover such universality relations between appropriately renormalised correlation functions in the 
continuum limit. However, note that with Wilson quarks it is a nontrivial matter to implement the 
chirally rotated boundary conditions ( Ol ), as it requires the fine tuning of a dimension 3 boundary 
counterterm (s. below). In this contribution we would like to address three questions: first, how 
difficult is it to implement the boundary conditions (|TTT|)? Second, given such an implementation, 



can we confirm relations following from universality, Eq. (1.3)? And finally, is automatic 0(a) 
improvement indeed realised? 

2. Lattice set-up 

The lattice set-up is taken from [g], with the fermion part of the action, 

T 

s /= (,4 LL#)(^+^+ ffl oW4 (2.i) 

jco=0 x 
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and the Wilson-Dirac operator: 

-Uo(x)P-\l/(x+aO) + (K+iy 5 z 3 P-)\i/(x), ifx = 0, 



-U Q (x)P-\i/(x + a6)+K\i/(x)-Uo(x-abyP + \l/(x-a6), ifO<x <T, (2.2) 
(K + iy 5 T 3 P+ )Y(x)- U (x - a6yP + y(x - a6) , for x = T . 



The time diagonal operator K is defined by 

K ¥ (x)= (l + ^{ a (V k + Vl) Yk - a 2 VtV k }) ¥ (x)+c m ja £ a^ v (x) V (x). (23) 

Finally, the counterterm action is specified by 

8&wy(x) = {S Xofl + K,T) (z f -l) + (d s -\)aD s y{x), (2.4) 

where the operator D 5 should reduce to Y%=\ Yk^k m the continuum limit [g]. With our conventions 
the tree-level coefficients are given by z*p = 1 and di ^ = 1 /2. While we set d s = <ij ' throughout, 
the finite renormalisation constant Zf must be determined non-perturbatively. 

3. Definition of correlation functions 

We need correlation functions for both the standard and the chirally rotated SF. In the standard 
SF we follow the conventions used in the literature jl0| ] by defining 

/x(*o) = -~ {x*H*)4 2h ) {p±) , M*>) = -\t {y( X HxWi 2h ) {p±) ■ (3-D 

Here the fields X" and Y k a stand for the quark bilinear fields, 

X=A ,V ,S,P, Y k =A k ,V k ,T k0 ,T kQ . (3.2) 

which are defined as usual, e.g. A^ 2 = YfiYsYfa- We also use the boundary-to-boundary corre- 
lators, 

For the chirally rotated SF we define correlation functions in the same way. Since the boundary 
conditions distinguish up and down type flavours we keep track of the flavour assignments by a 
superscript to the correlation functions. In order to avoid diagrams with disconnected fermion 
lines, we imagine a setup with 4 flavours, such that there are 2 up-type flavours and 2 down-type 
flavours. This greatly increases the flexibility when performing a chiral rotation, which can either 
rotate two flavours of the same type or two flavours of different types, while avoiding any Wick 
contractions corresponding to diagrams with disconnected fermion lines. The correlation functions 
in the chirally rotated set-up are denoted by gx and ly, as well as gi and l\. Their definition is such 



that universality implies the following relations from Eq. (1.3): 



f A = g »/ = -ig?f, f p = i gf=gf, k y=l^ = -ilf, (3.4) 
k T = ilf=Hf, f l=g f=gf, kl =lf = lf. (3.5) 
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T = L = 8, 13 = 6.0219, c sw = NP 
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Figure 1: Left panel: varying zj does not change the slope of g "f* (xq), which is proportional to hjpcac- Right 
panel: different tuning conditions for z/ yield O(a) differences. 



Here the flavour indices correspond to the quark bilinear operator being inserted, and the conven- 
tions for the quark boundary fields are taken from ref. [g]. All remaining correlation functions, such 
as /v or g{f' are expected to vanish by parity and flavour symmetries (as defined in the standard SF 
basis). 



4. Numerical simulations and results 

We have carried out a quenched simulation measuring the correlation functions for both the 
standard SF and chirally rotated SF on the same gauge configurations with vanishing gauge bound- 
ary fields. Both unimproved and non-perturbatively O(a) improved Wilson quarks were used [|TT|]. 
The chosen lattice sizes were (L/a) 4 with L/a = 8, 12, 16,24,32. The physical size of the lattice 



was kept fixed in terms of Sommer's scale ro [12], L/ro = 1.436, using the parameterisation [13] 



and the interpolation [|14|]. The critical mass m cr was set by requiring the PC AC mass to vanish in 
the standard SF. The simulation code is a customized version of M. Luscher's DDHMC code [0], 
and the simulations were run on PC clusters and on a BlueGene/L system. 

4.1 Tuning of z f, boundary conditions 

The tuning of z/ can be performed by requiring any ysT^odd quantity to vanish [g]. Ex- 
amples for such quantities are gp u (xq), or g^(xo). Fortunately, the tuning of the parameters mo 
and Zf is straightforward, as the respective tuning conditions are almost independent of each other 



(cf. also []16[]). Given the critical value m cr of mo, we have also checked that the differences Azf 
between z/ values obtained from different conditions vanish with a rate oc a, as expected (cf. fig. [[]). 
Given Z/(go) one then expects that the boundary conditions are correctly implemented up to 
cutoff effects. To test this hypothesis we have reverted the projectors Q± — > Q T in the bound- 
ary sources and indicate this change by a subscript "— " to the correlation functions. In the left 
panel of fig. |2| one indeed observes that the effect is very small and decreases towards zero. The 
corresponding result for the standard SF is very similar and given in the right panel of fig. ^[ 
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Figure 2: Check of boundary conditions, cf. main text for an explanation. 
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Figure 3: Universality checks: we consider ratios which should approach unity in the continuum limit. 



4.2 Universality relations 



The universality relations (3.4), (33b are expected to hold between renormalized correlation 



functions. For instance one should then find that the ratios 

-l 



gf{T/2)/^i) x A(r/2)/V/i 



/ P (r/2)/V7i 



-i 



(4.1) 



approach unity in the continuum limit. As seen in fig. ||, this is well satisfied within errors. 

4.3 Automatic 0(a) improvement 

Automatic bulk O(a) improvement relies on the fact that the bulk O(a) effects are located in 
75 T 1 -odd correlation functions. By projecting on the 75 T 1 -even correlators one thus gets rid of 
the bulk O(a) effects. Here we study the 75 T 1 -odd correlators and verify that these vanish in the 
continuum limit with a rate oc a, as can be seen in fig. ^| for 1%" (T/2). The corresponding standard 
SF correlator is also shown and vanishes exactly after gauge average, as expected due to parity 
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Figure 4: Check of automatic O(a) improvement, cf. main text for details. 



and flavour symmetries. Another example is the counterterm contribution oc c\(a/L) needed to 
improve correlation functions of the axial current [0]. The result in this case is shown in the right 
panel of fig. ||. While the continuum limit vanishes in the chirally rotated SF, it is finite in the 
standard SF. Hence its contribution to axial current correlators is of 0(a 2 ) and 0(a), respectively, 
thereby confirming the expectation from automatic O(a) improvement. 



4.4 Determination of Z, 



A.V 



Having checked universality we may turn the tables and use universality to determine a number 
of finite renormalization constants which are usually determined from chiral Ward identities. For 
instance, the continuum relations ( fO| ) imply that Za,v can be determined by the ratios, 

Z A = igtf(T/2)/fl?(T/2), Zy = #(T/2)/itf(T/2), (4.2) 

where denotes the conserved vector current. The results in fig. || are subject to an 0(a 2 ) uncer- 



tainty, which perfectly explains the discrepancy with the Ward identity results of ref. [17]. 



5. Conclusions 

We have presented a successful implementation of the chirally rotated Schodinger functional. 
Universality could be confirmed by comparing with standard SF correlation functions and auto- 
matic O(a) improvement has been verified. In the future, we expect this framework to yield better 
controlled continuum extrapolations for step-scaling functions. Furthermore it provides new meth- 
ods to determine finite renormalizaiton constants and O(a) improvement coefficients for gauge 
theories with Wilson- type fermions. 
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